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Molecular Hydrogen in a Damped Lyman-a System at Zabs = 4.224 

C. Ledoux^, P. Petitjean^'^, R. Srianand^ 

O . ABSTRACT 
O ■ 
(N ; 

' We present the direct detection of molecular hydrogen at the highest red- 

pL( i shift known today (2;abs = 4.224) in a Damped Lyman-a (DLA) system toward 

Q\ I the quasar PSS J 1443+2724. This absorber is remarkable for having one of 

^ I the highest metallicities amongst DLA systems at Za^s > 3, with a measured 

>• I iron abundance relative to Solar of —1.12 ± 0.10. We provide for the first time 

^ I in this system accurate measurements of Nl, Mgll, Sll and An column den- 

' sities. The sulfur and nitrogen abundances relative to Solar, —0.63 ± 0.10 and 

Q I —1.38 ± 0.10 respectively, correspond exactly to the primary nitrogen produc- 

^ ' tion plateau. H2 absorption lines are detected in four different rotational levels 

^ ' (J = 0, 1, 2 and 3) of the vibrational ground-state in three velocity components 

Ph' with total column densities of logA^(H2) = 17.67, 17.97, 17.48 and 17.26 respec- 

O ' tively. The J = 4 level is tentatively detected in the strongest component with 

' logA^(H2) ~ 14. The mean molecular fraction is log/ = —2.38 ± 0.13, with 

; / = 2A^(H2)/(2Ar(H2) + A^(Hl)). We also measure logiV(HD)/iV(H2) < -4.2. 

^ ^ ■ The excitation temperatures Tqi for the two main components of the system are 

^ ' 96 and 136 K respectively. We argue that the absorbing galaxy, whose star- 

ch ■ formation activity must have started at least 2 — 5 x 10^ yrs before z = 4.224, is 

in a quiescent state at the time of observation. The density of the gas is small, 
nn < 50 cm~^, and the temperature is of the order of T ~ 90 — 180 K. The 
high excitation of neutral carbon in one of the components can be explained if 
the temperature of the Cosmic Microwave Background Radiation has the value 
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expected at the absorber redshift, T — 14.2 K. These observations demonstrate 
the feasibihty to study H2 at the highest redshifts provided high enough spectral 
resolution and good S/N ratio are achieved. 

Subject headings: cosmology: observations - galaxies: formation - ISM: molecules 
- quasars: absorption lines 

1. Introduction 

Although molecular hydrogen, H2, is the most abundant molecule in the Universe, it 
is very difficult to detect directly. Emission lines from the tracer molecule CO are usually 
used instead but, at z > A, only the brightest objects can be detected with current facilities 
(e.g., Omont et al. 1996; Solomon & Vanden Bout 2005). For the time being, the only 
way to detect H2 directly at high redshift is to search for the absorption signature of the 
molecular UV transition lines from Damped Lyman-a (DLA) systems observed in the spectra 
of quasars or those of 7-ray burst afterglows. 

DLA systems are characterized by very large neutral hydrogen column densities, iV(H l) > 
10^° cm~^, that are typical of sighthnes through local spiral galaxies. Inferred metallicities 
typically vary between [Zn,S/H] = —2 and -0.5 at 2 < Zg^s ^ 3 and are smaller at larger 
redshifts (e.g., Prochaska et al. 2003). The dust content is less than or about 10% of what 
is observed in the Galactic ISM for similar Hi column densities (Ledoux et al. 2003). The 
detection of H2 molecules in high-redshift DLA systems through Lyman- Werner band ab- 
sorption is observationally challenging due to the presence of the Lya forest. High spectral 
resolution, good S/N ratio and large wavelength coverage arc simultaneously needed to de- 
tect the H2 lines and avoid spurious detections (e.g., Lcvshakov & Varshalovich 1985; Ge & 
Bcchtold 1999; Petitjean et al. 2000; Levshakov et al. 2002; Ledoux et al. 2002, 2003; Gui 
et al. 2005). 

We are conducting a survey for molecular hydrogen in DLA systems at high redshift 
(^^abs > 1-8; see Ledoux et al. 2003), using the Ultraviolet and Visible Echelle Spectrograph 
(UVES; Dekker et al. 2000) at the European Southern Observatory's Very Large Telescope, 
down to a detection hmit of typically A'"(H2) = 2 x 10-*^^ cm~^. Out of the 75 systems 
observed up to now, 14 have firm detections of associated H2 absorption lines. In this 
Letter, we present the highest redshift detection in our sample, in the z^^s — 4.224 DLA 
system toward the quasar PSS J 1443-1-2724. We describe the observations in Sect. 2 and 
the analysis of the data in Sect. 3. We conclude in Sect. 4. 
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2. Observations 

The quasar PSS J 1443+2724 {zem = 4.42, r ~ 19.3) was observed with UVES in service 
mode on March 2-3 and 17-19, 2004. The red spectroscopic arm of UVES was used in the 
standard configuration with the central wavelength adjusted to 580 nm. The wavelength 
coverage was 478—681 nm with only a small gap, between 575 and 584 nm, corresponding 
to the physical gap between the two red arm CCDs. The CCD pixels were binned 2x2 
and the sht width was fixed to 1", yielding, under the typical 0'.'7 seeing conditions achieved 
during the observations, a resolving power R 48000. The total on-source integration 
time was 7.3 h, split into five exposures. The data were reduced using the UVES pipeline 
(Ballester et al. 2000) which is available as a dedicated context of the ESO MIDAS data 
reduction system. The wavelength scale of the spectra reduced by the pipeline was converted 
to vacuum-heliocentric values and individual exposures were scaled, weighted and combined 
altogether to produce the final object fiux and variance spectra. 

3. Analysis 

3.1. The high metalUcity and the star- formation history 

The metaUicity of the DLA system toward PSS J 1443-1-2724 is the highest among the 
ten z > A DLA systems investigated up to now (Prochaska et al. 2003). Total integrated 
column densities and metaUicities for Nl, Mgll, Sll, An and Fell, obtained from Voigt- 
profile fitting of the observed transition lines, arc given in Table 1. With the exception of 
Fe II, these ions are observed here for the first time. We derive metaUicities relative to Solar 
of [S/H] = -0.63 ± 0.10 and [Fe/H] = -1.12 ± 0.10. The latter value is consistent with the 
previous measurement by Prochaska et al. (2001) but we note that the newly determined 
Hi and Fell column densities are both larger by about 0.15 dex. Even though the depletion 
factor is moderate (i.e., [S/Fe] = -1-0.49), because of the high metaUicity the dust-to-gas ratio 
measured in this system (logK = —0.8) is similar to that of other systems in which H2 is 
detected (see Fig. 15 of Ledoux et al. 2003). Substantial star-formation activity must have 
taken place in the past history of this DLA system. Note that the age of the Universe at 
z = 4.22 is 1.5 Gyr (for Qa = 0.73, = 0.27 and Hq = 71 km s-^ Mpc"^). 

Nitrogen is thought to have both primary and secondary origins depending on whether 
the seed carbon and oxygen nuclei are produced by the star itself or are already present 
in the ISM from which the star formed. In the case of secondary nitrogen production, the 
ratio of nitrogen to oxygen abundances increases with increasing oxygen abundance, whereas 
for primary production this ratio remains constant and nitrogen tracks oxygen (this is the 
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so-called primary plateau). In the low-metallicity Hll regions of nearby galaxies, [N/0] 
abundance ratios tend to lie around the primary plateau (with [N/0] = —0.73 ± 0.13) 
for [0/H] < —1, and have a secondary behaviour for [0/H] > —1 (van Zee et al. 1998; 
Izotov & Thuan 1999). Similarly, most DLA systems show a plateau at [N/0] ~ —0.9 
(see Centurion et al. 2003; Pettini et al. 2002). The abundance of nitrogen relative to 
sulfur at z^-bs = 4.224 toward PSS J 1443+2724, [N/S] = -0.75, corresponds exactly to the 
primary nitrogen production plateau (assuming [0/S] = 0). If we assume the widely accepted 
hypothesis that primary nitrogen is produced in intermediate-mass stars (4 < Mq < 8) 
during the Asymptotic Giant Branch phase (Henry et al. 2000), then star- formation activity 
must have started in this system at least 2 — 5 x 10^ yrs before z — 4.224 accounting for the 
characteristic lag time for intermediate-mass stars to eject nitrogen. This is similar to what 
is observed in the Zabs = 4.383 system toward BRI 1202-0725 (D'Odorico & Molaro 2004). 

Argon can be ionised in H I regions and is therefore a tracer of the radiation field. In 
the presently studied system, we find [Ar/S] = —0.30, which is higher than in most other 
DLA systems (Vladilo et al. 2003). This is consistent with the associated galaxy being in 
a quiescent state at the time of observation, as also inferred below from the analysis of H2 
lines. 



3.2. Moleculcir hydrogen 

The H2 transitions are numerous and the spectral resolution of our data is high enough 
to allow unambiguous detection and detailed analysis of H2 (see Fig. 1). Great care has 
been exercised when fitting the absorption lines with multi-component Voigt profiles. H2 is 
detected in four different rotational levels (J = 0, 1,2 and 3) of the vibrational ground-state 
in three velocity components with total column densities of log A^(H2) = 17.67, 17.97, 17.48 
and 17.26 respectively (see Tables 1 and 2). The J = 4 level is tentatively detected in the 
strongest component with logA'^(H2) ~ 14. Blending is heavy however and this could be 
an upper limit. The detection of H2 confirms the previous claim that most of the gas in 
this DLA system is cold (Howk et al. 2005). The total H2 column density is one of the 
largest observed in DLA systems and the molecular fraction is log / = —2.38 ± 0.13, with 
/ = 2A^(H2)/(2A^(H2) + iV(Hl)). HD is not detected down to 3a limits of logiV(HD) < 13.6 
and 13.7 for J = and 1 respectively, implying that [HD/H2] < —4.2. 

The excitation temperatures Tqi for the two main components of the system are 96 and 
136 K respectively (see Table 2). The excitation temperatures for higher rotational levels 
are slightly larger, possibly as a consequence of the levels being populated by radiative or 



- 5 - 



formation pumping. By equating H2 formation and destruction rates, we write: 

SIn{R2) = Rn{R)n (1) 

where n = n(H)+2n(H2) is the particle density, / is the H2 dissociation rate in the unshielded 
UV background (in s^^) and S = [A^(H2)/10^^ cm~^]~°'^^ is the correction factor for shielding 
(Draine & Bertoldi 1996; Hirashita & Ferrara 2005). i? = 4.1 x lO-^^fT/lO^ Kf^{Z/ZQ)D 
cm^ s~^ is the H2 formation rate per unit volume and time, Z/Zq = 0.23 is the metallicity 
relative to Solar and D = 0.7 is the fraction of metals into dust. 

The H2 J = 4 level is mostly populated by pumping of lower rotational levels by ab- 
sorption of UV photons from the ambient radiation field and subsequent radiative cascade, 
and by direct formation of the molecule in this excited state. Following the simple procedure 
introduced by Jura (1975), we write the equation describing the equilibrium of the J — A 
level population as: 

9.1p4,0'5/n(H2, J = 0) + 0.19i?n(H)n = ^(4 ^ 2)n(H2, J = 4) (2) 

where p4 = 0.26 is the pumping efficiency into the J = 4 level from the J = level and 
A{A 2) = 2.8 X 10""^ s~^ denotes the spontaneous transition probability from J = 4 to 
J = 2. 

Combining the above equations and using the total column densities given in Table 1, 
we find Rn = 4 x 10~^^ s~^. This is probably an upper limit given the possibifity that 
A^(J = 4) be an upper limit. Assuming a temperature of T" = 150 K, we derive a density of 
TT-H < 50 cm~^. This value could be slightly larger in case of any a-element enhancement. The 
corresponding H2 dissociation rate is / < 3 x lO^^*' s~^. This is of the order of or smaller 
than the highly variable photo-dissociation rate observed in the ISM of the Galaxy (e.g., 
Hirashita & Ferrara 2005). This probably indicates, as expected above, that there is little 
on-going star formation in the galaxy associated with this DLA system. This is consistent 
with the fact that no emitting counterpart has been found using the Lyman-break techniques 
down to L ~ (Prochaska et al. 2002) or in Lyman-a emission down to 2 x 10"^^ erg 

cm~^ s~^ from 2D-spectroscopy at the Canada- Prance-Hawan Telescope (Ledoux 1999). 



3.3. The Cosmic Microwave Background temperature 

Absorption fines from the ground-state and first fine-structure excited level of C I are 
detected (see Fig. 2) at redshifts of Zabs — 4.22379 and 4.22413, sfightly larger than those of 
the molecular components. The Doppler parameters are larger as well indicating that the C I 
gas spans a slightly larger velocity range than the H2 gas (see also Ledoux et al. 2003). We 
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measure [C I*/ C l] = — 0.09±0.05 and — 0.21±0.04 in each of the two components respectively. 

This is surprisingly small at high redshift. Indeed, the value expected for excitation by the 
Cosmic Microwave Background Radiation (CMBR) field alone, with a temperature of 14.2 K 
at z = 4.22, is —0.15 (Silva & Viegas 2002). In addition, part of the excitation can be due to 
UV or IR pumping although this is probably negligible here. In Fig. 3, we plot the expected 
values for [Cl*/Cl] as a function of temperature of the cold gas. The first Cl component 
has a density (nn < 25 cm~^) that is similar to what is measured in the H2 gas and is 
thus probably predominantly associated with it. The second C I component probably has a 
lower density and the corresponding [C I*/ C l] ratio is consistent with what is expected for 
excitation by the CMBR field alone. 

4. Conclusion 

The direct detection of H2 at high redshift is a unique tool for studying the ISM in 
the remote Universe. We have detected molecular hydrogen in one of the highest redshift 
DLA systems known to date (zabs = 4.224). We have shown that the gas is probably the 
left-over of intense star-formation activity that happened at least 2 x 10^ yrs before the time 
of observation and that the associated object is in a quiescent phase. The UV radiation field 
is of the order of or smaller than that in the Galaxy and the gas is cold (T ~ 90 — 180 K). The 
[Cl*/Cl] ratio in one of the components of the system is consistent with what is expected 
for excitation by the CMBR field at this redshift (T = 14.2 K). This study demonstrates the 
importance of this kind of observations for our understanding of how star formation proceeds 
at high redshift. The prospect of using 7-ray burst afterglows as background sources for 
similar studies is exciting and will boost this field in a near future (Draine & Hao 2002). 

This work is based on data collected under prog. 072.A-0346 of the European Southern 
Observatory with UVES mounted on the 8.2 m Kueyen telescope at the Paranal Observa- 
tory, Chile. We thank the referee, Paolo Molaro, for fruitful comments. PP and RS acknowl- 
edge support from the Indo-Prench Centre for the Promotion of Advanced Research (Centre 
Franco- Indien pour la Promotion de la Recherche Avancee) under contract No. 3004-3. 
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Table 1: Total column densities of different species in the DLA system at Zahs = 4.224 toward 
PSS J 1443+2724 
Species Transition log N ± criog at Metallicity 

lines relative to Solar ^ 



Hi 


1025,1215 


20.95 ±0.10 




H2 


J=0 lines 


17.67 ±0.06 




H2 


J=l lines 


17.97 ±0.05 




H2 


J=2 lines 


17.48 ±0.12 




H2 


J=3 lines 


17.26 ±0.17 




H2 


J=4 lines 


14.02 ±0.14 




Ci 


1277,1280 


13.12 ±0.02 




Ci* 


1277,1279 


12.97 ±0.03 




Ci** 


1277,1280 


< 12.5 




Nl 


953.4,1134.1 


15.52 ±0.01 


-1.38±0.10 


Mgll 


1239,1240 


15.98 ±0.01 


-0.55±0.10 


Sii 


1250 


15.52 ±0.01 


-0.63±0.10 


An 


1066 


14.42 ±0.01 


-0.93±0.10 


Fell 


1081 


15.33 ±0.03 


-1.12±0.10 



^ Solar abundances adopted from Grevesse & Sauval (2002). 

Note: errors in the column densities of metals correspond to the rms errors from fitting the 
Voigt profiles. They do not include the uncertainties related to the continuum placement. 
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Table 2: Column densities and excitation temperatures from Voigt- profile fitting of H2 molec- 
ular lines 



-^■abs 


Level 


Ir-vrr AT 

log iV 

iCTlogTV 


± CTf, 

[km s^"^] 


ex 

[K] 


T 

"ex 


4.22371 


J=0 


1741+0.02 

' -^^^-0.05 


0.8 ±0.3 








J=l 


1 7 f:q+0.03 


5) 


90-107 


0-1 




J=2 


i7.i5in^ 


5) 


200-259 


0-2 




J=3 


16.53l°il 


)) 


156-302 


1-3 




J=4 


< 13.6 ^ 


)) 






4.22401 


J=0 


7 04+0.09 


1.1 ±0.3 








J=l 


1 7 7c;+0-06 


)) 


104-179 


0-1 




J=2 


17 2Ci+0.13 




218-377 


0-2 




J=3 


17 91+0.13 
^ ' •^-'--0.13 


)) 


337-530 


1-3 




J=4 


14.UZ_o.o7 


)) 






4.22416 


J=0 
J=l 
J=2 
J=3 
J=4 


14.34+n^ 

14.86lH^ 

14.6llH^ 
14.181H^ 

< 13.6 ^ 


1.7 ±0.5 

5) 
5) 
)) 
?) 







^ 3(7 detection limit. 

Note: errors in the column densities correspond to the range 6 ± (76 of Doppler parameters 
and not to the rms errors from fitting the Voigt profiles. 
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Fig. 1. — Velocity profiles of selected transition lines from the J = 0, 1, 2, 3 and 4 rotational 
levels of the vibrational ground-state Lyman and Werner bands of H2 at z^hs = 4.224 toward 

Tiriri T -1 A A I n'-7n A ml 1_ ^ T 7"^ • _-J_ J21 ^ • ^ ^1 J_l 1 ^1 ^ • J_l_ 
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Fig. 2. — Velocity profiles of Cl and Cl* transition lines at ^abs — 

4.22379 and 4.22413 

toward PSS J 1443+2724. The best Voigt-profile fit is superposed on the observed spectrum 
with vertical dotted lines marking the location of the two detected components of the system 
and a third, possible component. 
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Fig. 3. — [Cl*/Cl] ratio versus temperature. Filled circles with error bars are the mea- 
surements in the ^abs = 4.22379 and 4.22413 components toward PSS J 1443+2724. LTE 
predictions for two different values of the particle density are shown by solid lines. The value 
expected for excitation by the CMBR field alone is shown by the horizontal long-dashed line. 



